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Abstract 

Time dependent entropy of harmonic oscillator with time dependent mass and frequency are 
investigated. The joint entropy so called Leipnik's entropy is calculated by using time dependent 
wave function obtained by the Feynman path integral method. It is shown that, Leipnik's entropy 
fluctuates with time. However in constant mass and time dependent frequency case, entropy 
increases monotonically with time. 
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I. INTRODUCTION 



The information entropy plays a major role in a stronger formulation of the uncertainty 
relations [1]. This relation may be mathematically defined by using the Boltzmann-Shannon 
information entropy and the von Neumann entropy. For both open and closed quantum 
systems, different information-theoretical entropy measures have been discussed |2|, y, ]4\. 
The joint entropy [5, f| can also be used to explain the properties of the loss of information 
evolving pure quantum states & The joint entropy of the physical systents were conjectured 
by Dunkel and Trigger [8] in which their systems named MACS (maximal classical states). 
The Leipnik entropy of the simple harmonic oscillator was determined not monotonically 
increase with time (9|. In this work, we give a uniform description of the complete joint 
entropy information for time dependent entropy of harmonic oscillator with time dependent 
mass and frequency. The study of harmonic oscillators with time dependent frequencies or 
with time dependent masses (or both simultaneously) has attracted considerable interest in 
past few years but the investigation of joint entropy of this system have not been studied 
yet. The time-dependent harmonic oscillator has invoked much attention because of its many 



app lications in different areas of physics, such as quantum optics and plasma physics [10, 
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This paper is organized as follows. In section II, we explain fundamental definitions 
needed for the calculation. In section III, we get the results for harmonic oscillator with 
time dependent mass and frequency. Moreover, we obtain the analytical solution of kernel 
and using this, wave function in both coordinate and momentum space and its joint entropy 
were calculated. In this section, we also investigated same quantities for harmonic oscillator 
with strongly pulsating mass and inverse square time dependent frequency. Finally, we 
present the conclusion in section IV. 

II. FUNDAMENTAL DEFINITIONS 

We consider a classical system with d = sN degrees of freedom, where N is the par- 
ticle number and s is number of spatial dimensions 8|. The density function g(x,p,t) = 
g(xi, ...,Xd,Pi, ...,Pd,t) which is non-negative, time dependent phase space density function 



2 



of the system is assuming to be normalized to unity, 



dxdpg(x,p, t) — 1. 



(1) 



The Gibbs-Shannon entropy is described by 

1 



(2) 



where h = 2irh is the Planck constant. Schrodinger wave equation with the Born interpre- 



tation 



141 ] is given by 



dip 

in—— = Hib. 
dt Y 



(3) 



The quantum probability densities are defined in position and momentum spaces as \ip(x, t)\ 2 
and \ip{p, t)\ 2 , where \ip(p,t)\ 2 is given as 

dxe~ ipx/h 



Leipnik proposed the product function as 



{2iih) d i'- 



-i>(x,t). 



(4) 



gj(x,p,t) = \tP(x,t)\ 2 \tP(p,t)\ 2 > 0. 



(5) 



Substituting Eq. (5) into Eq. (2), we get the joint entropy Sj(t) for the pure state ip(x,t) 
or equivalently it can be written in the following form 



S&) = - J dxmx,t)\ 2 \nmx,t)\ 2 - I dp\^{p,t)\Hn\^{p,t)\ 2 
- \nh d . 



(6) 



We find^time dependent wave function by means of the Feynman path integral which has 
form 



K(x",t";x',t" 



x"=x(t") 



Dx{t)e 



x'=x(t') 

Dx(t)e^^' L[x ' ± ' t]dt . 



(?) 



The Feynman kernel can be related to the time dependent Schrodinger's wave function 

oo 

K(x", t"; x', t') = ]T O- (8) 



n=0 
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The propagator in semiclassical approximation reads 

i d 2 



K(x",t";x',t') 



Scl{% jt j X , t ^ 



1/2 jrs cl (x",t"y,t') 



(9) 



.2nh dx'dx" 

The prefactor is often referred to as the Van Vleck-Pauli-Morette determinant [l6j, [l7j. The 
F(x", t"; x', t') is given by 

i d 2 i J ^ 



F(x"X;x',t') 



2nTi dx'dx 



j— Sd(x , t ; x , t ) 



(10) 



III. HARMONIC OSCILLATOR WITH TIME DEPENDENT MASS AND FRE- 
QUENCY 

The Lagrangian of the harmonic oscillator with time- dependent frequency and mass are 
given by 

L = ^m(t)x 2 - ^m(t)u 2 (t)x 2 , (11) 

where u>(t) and m{t) are, respectively, the frequency and mass associated with oscillator, 
and which are arbitrary real functions of time. The classical equation of motion for the 
Lagrangian is 

x + 2^x + u 2 (t)x = 0, (12) 
V 

where r](t) = \fm. The solution to the equation of motion is given by 

a(t) 



x(t) 



T)(t) 



A cos 7(i) + B sin j(t) 



(13) 



where = p(t), ^(t) refers to the amplitude and phase of classical oscillators and A and 
B are constants. The function p(t) have to satisfy the following equation: 



p+-p + U j 2 (t)p(t) 



m m 2 p 3 



(14) 



The constant A and B in Eq.(13) can be determined by using the boundary conditions of 
x{t') = x' and x(t") = x" . The A and B coefficients yield as 

1 



and 



.4 



B 



sin (7" — 7') \p 



x 



X 



— sin 7 sin 7 

p" 



x 



X 



, - cos 7 cos 7 

sin (7" — 7') V p" p' 



(15) 



(16) 
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Substituting A and B into Eq.(13), the classical path that connects the point of (x',t') and 
(x",t") can be written as 

|_^(^3in(y- 7 )-^sin(y'- T )). (17) 



d sin(7" — 7') V p' 



Substituting the classical paths into action function, the classical action becomes 

1 1 1 I 111 III / 111 \ I I 

1 " ' 2 \p'/ 2 \p"* p" 



&)-^6Hfe + ^H^>- 



x"x' 



p"p 



2 \p 
7 csc(t"-t'). 



By substituting the above action into Eq.(lO), the factor can be obtained as 



F(x",t";x',t') 



1 



-,1/2 



(18) 



(19) 



-2iixTip' p" sin(7" — 7' 

Therefore the propagator for the harmonic oscillator with a time dependent mass and fre- 
quency can be obtained by 

1 



K(x",t";x',t') 



2mhp' ' p" sin(7" — 7') 
exp 



1 1/2 


" i ' 


exp 


I2h[ 



m"x" 2 



10 • I 

m x ( p 



2 \p" 



x cos(7" - 7') 



2^sin(7"-Y) \p" 2 p' 2 
2x"x' 



P'P" 



(20) 



By the use of the Mehler-formula 

-C- 



e-^y 2 J2^) 2 H n (x)HM 



:CXp| 



n=0 



VT- 

(x 2 + y 2 )(l + z 2 ). 



4xyz 
2{l-z 2 ) 



2(1 -z 2 ) J ' 

where H n is Hermite polynomials, we can write the Feynman kernel in form 



(21) 



K(x", t"; x', t') 



2 n n\ 



1 



to • t 

mx / p 



.2%TxUp' p" sin(7" — 7' 
exp 



1/2 



exp 



% \m"x" 2 /p"\ 



x H„ 



2 \p> 

V /|^]exp(-,7(n+^)). 



1 /x" 2 x' 2 

~ 2h \~jp* + ^ 



2hl 2 \p" 



x 



(22) 



Comparing the kernel in Eq. (22) with Eq. (8), the wave function for harmonic oscillator 
with time dependent mass and frequency can be found as 



V n (x,t) = exp[ia n (t)] 



1 



x H„ 



^/ 2 n 1/2 n\2 n p 
1\ 1/2 x 



1/2 



exp 



im(t) / p 



+ 



2h \p m(t)p 



x 



x 



\h) ~p. 



(23) 
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where the phase functions a n (t) are described by 

a n (t) = -(n + 7i) [ u>\ 2 dt '- 
\ 2/ J m(t')p 2 



(24) 



This result is also briefly obtained before by using the Feynman Path Integral method in 
Ref.P 

dependent wave function of ground state is 



ggral 

. It is also in agreement with ones obtained in Refs. [lsj and [h]]. The time 



Vo(x,t) = exp[ia (t)} 



I 1 


1/2 r 




exp 



im(i) / p 



+ 



x 



(25) 



2h \p m(t)p 2 , 

Note that when m(t) — > m, oj(i) — > co> and p(t) — > p = const = 1/y/mujQ, the solution in 
Eq.(23) becomes the solution for the time independent harmonic oscillator of mass m and 
frequency u . 



A. Time dependent mass and constant frequency 



In this calculation, we choose harmonic oscillator constant frequency and time dependent 
mass which is called pulsating mass. For instance, it has been shown that the Lagrangian 
describing the problem of a Fabry- Perot cavity in contact with a heat reservoir assumes the 
form of constant frequency and time dependent mass 20] . The time dependence of the mass 
can be taken as m(t) = mcos 2 (z/t) where v is described as frequency of pulsating mass. The 
Lagrangian of the system is 



L = -m cos 2 (ut)x 2 — -mcos 2 (i/f)wV. 



(26) 



The kernel for a harmonic oscillator with strongly pulsating mass can be obtained from Eq. 
(20) as 



K{x\t"-x\t') 



mQ cos vt' cos vt" 



2mhsinVL(t" - t') 
cos 2 vt' tan vt! x 12 



1/2 



exp 



imv 



exp 



2h 
imQ 



27tsin(fit" -f) 



2 4. i // 112 

cos vtanvt x — 



2 i // 112 i 

cos vt X + 



+ cos vt'x' 2 cos Q(t" - t') - 2 cos vt' cos vt"x"x' 



(27) 

where we use Vt 2 = u 2 + v 2 . Using Eq. (21) Mehler formula, and then the time dependent 
wave function was obtained as 



x K 



m(t)n 



1/4 



7r7m! 2 2 2ri J 
m(t)tty/ 2 



exp 



im{t) 
L 2h 



z/tan vt + iQ ) x 



x 



x 



(28) 



where the phase functions (3 n {t) is described as 

Pn(t) = (n + l/2)Qt. 



(29) 



Note that when v — > and m(t) — > m = const, the solution in Eq.(28) reduces to the 
solution for the time independent harmonic oscillator of mass m and frequency Uq. The 
ground state wave function yields as 

~im(t) 



im(t)Q- 


1/4 r 




exp 



V (x,t) = exp[-i(l/2)0^ ^ j exp ( u tan ut + iVt ) x 2 . (30) 

The density of probability in coordinate space is 

IT . ,, 2 rm(t)0-|V2 r m(t)VL 2 i 

*o M 2 = -H- ex p ■ 31 

L 7171 J L ft J 

In Fig{TJ the density of probability in coordinate space is shown. The density of probability 
in momentum space can be easily calculated as 



l*o(p,t)| 5 



n 



fmrniy 2 tan z/t + fi 2 ) 



1/2 



exp 



p 2 ft 



mh(is 2 tan z/i + Q 2 ) 



(32) 



The joint entropy for ground state from Eq. (6) becomes 



1 

2 L 



In 



V4/ Wfi,^ 



+ A/m(t)fi( 1 - In 



(33) 



7i7i(z/ 2 tan 2 z/t + fi 2 ). 

The joint entropy of harmonic oscillator with time dependent mass and constant frequency 
is plotted in FigfSJ The small v and large v cases are shown in Figj3] and FigJH respectively. 



B. Time dependent frequency and constant mass 

The Lagrangian of the harmonic oscillator with time-dependent frequency is given by 

(34) 

The time dependent frequency is described by 

u/(t) = (35) 



t -2 1 /j.\2 2 

L = —mx mult) x . 

2 2 w 
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In this case, we denned the following quantities 



Pit) 



y/mu}(t) 



and 



lit) = - — . 



(36) 



(37) 



By substituting Eqs. (35), (36) and (37) into the general kernel Eq. (20), the kernel for a 
harmonic oscillator with the inverse square time dependent frequency can be derived as 



K(x",t";x',t') 



mujQ 



x exp 



2iirht't" sin uq [ ,,,,, 
imujQ 



t"-t' 



1/2 



exp 



//9 /9 

% /mx mx 



2h\ t" 



r 



X 



fx" 2 x' 2 \ ft" -t' 



-2h sin uj {^—± 
Using Eq. (21) Mehler formula, we obtain the kernel 

K(x",t";x',t') 



H -r I COSWo 



l\f' 2 t' 2 J 



t'f 



2x'x'' 
ft" 



' muj " 


1/2 r 


Inhff- 


exp 



It 1 ? /9 

% /mx mx 



x 



x exp 



muo 
~2~h 



JI2 



fx 

LV? 2 1 f 2 ) 



x' 2 Y N 



mujQ x 
~h~l" 



mojQ x 



exp(*f (n + 1/2)) 



Using the Eq. (8), the time dependent wave function was derived as 



tf„(x,t) 



1 



V2 n n\t 
x H 



mujQ 
7rh - 



1/4 



exp I I n + 



2/ t / 



exp 



2frA t J 



x 



mu a; 
^ 1 



The time dependent wave function for ground state is given by 



mujQ 
iM 2 



1/4 



exp 



\2t J 



The density of probability in coordinate space is 

|^o(x,t)| 2 



exp 



2ht\ t 



~muJo~ 


1/2 


mujQ 2 ' 


WntA 


exp 


[ w x \ 



(38) 



(39) 



(40) 



(41) 



(42) 

The density of probability for coordinate space is shown in Figj5j The density of probability 
in momentum space can be easily calculated as 



l*o(p,*)| 5 



u t 2 



.(u)q + t 2 )hnm. 



1/2 



exp 



p 2 u t 2 



(43) 
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Note that when t — > 1 and u(t) ujq = const, the solution of Eqs.(42) and (43) become 
the solution for the time independent harmonic oscillator of mass m and frequency ujq. The 
joint entropy for ground state from Eq. (6) becomes 



The contour and joint entropy of harmonic oscillator with inverse square time dependent 
frequency are shown in Figs. [6] and [7J respectively. It is important that Eq. (44) is in 
agreement with following general inequality for the joint entropy: 



originally derived by Leipnik for arbitrary one-dimensional one-particle wave functions^, |8j. 

IV. CONCLUSION 

We have investigated the joint entropy for explicit time dependent solution of one- 
dimensional harmonic oscillators with time dependent frequency and mass. The time de- 
pendent wave function is obtained by means of Feynman Path integral technique. In the 
time dependent strongly pulsating mass, harmonic oscillator with constant frequency case, 
we have found that the joint entropy fluctuated with time and frequency. It is seen that, 
the joint entropy has harmonic behavior. This result indicates that the information pe- 
riodically transfers between harmonic oscillators with strongly pulsating masses. On the 
other hand, inverse square time frequency case, the joint entropy of harmonic oscillator 
with time dependent frequency shows remarkable monotonically increase with time. In the 
u(t) — > u — constant and v — * 0, case theses results agree with time independent harmonic 
oscillator. It also depends on choice of initial frequency as seen in Fig. [7J 
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FIG. 1: \^o(x,t)\ 2 versus time and coordinate for strongly pulsating mass and constant frequency. 
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FIG. 2: The 3D graph of joint entropy of harmonic oscillator with time dependent mass and 
constant frequency. 
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FIG. 3: The joint entropy of harmonic oscillator with time dependent mass and constant frequency 
versus small v. 
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FIG. 4: The joint entropy of harmonic oscillator with time dependent mass and constant frequency 
versus large v. 
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FIG. 6: The counter graph of the joint entropy of Harmonic oscillator with inverse square time 
dependent frequency and constant mass. 



16 




FIG. 7: The 3D graph of the joint entropy of harmonic oscillator with inverse square time dependent 
frequency 
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